Energy gap opening in submonolayer lithium on graphene: Local density functional 

and tight-binding calculations 
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The adsorption of an alkali-metal submonolayer on graphene occupying every third hexagon of 
the honeycomb lattice in a commensurate arrangement induces an energy gap 

in the spectrum of graphene. To exemplify this type of band gap, we present ab initio density 
functional theory calculations of the electronic band structure of CeLi. An examination of the 
lattice geometry of the compound system shows the possibility that the nearest-neighbor hopping 
amplitudes have alternating values constructed in a Kekule-type structure. The band structure of 
the textured tight-binding model is calculated and shown to reproduce the expected band gap as 
well as other characteristic degeneracy removals in the spectrum of graphene induced by lithium 
adsorption. More generally we also deduce the possibility of energy gap opening in periodic metal 
on graphene compounds C X M if x is a multiple of 3. 

PACS numbers: 73.22.-f, 73.20.At, 81.05.Uw 
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I. INTRODUCTION 

The isolation of graphene, a honeycomb lattice of car- 
bon atoms, and the observation of the electric- field ef- 
fect in the nanostructured samples, deposited on oxidized 
silicon surface, have renewed interest in the electronic 
properties of this two-dimensional carbon material^ The 
breakthrough discovery made possible the experimen- 
tal observation of theoretically predicted exotic physics 
of graphene, including an anomalous half-integer quan- 
tum Hall effect, non-zero Berry's phase, and minimum 
conductivity^^ 

The exceptional properties of graphene, two- 
dimensional structure and room-temperature high 
mobility, make it an ideal material for carbon-based 
nanoelectronics envisaged in the future nanotechnolo- 
gies. Various applications depend on the possibility 
of fabricating a graphene-based field effect transistor 
(FET), including transistors to be used in nanoscale 
electronic circuits, light emitters and detectors, and 
ultra-sensitive chemical sensors and biosensors 4. Being 
a zero-gap semiconductor, however, graphene cannot 
be used directly as the conducting channel in FETs. 
Discovering how to generate band gaps is therefore 
a question of fundamental and applied significance. 
Several possibilities exist. Lateral confinement of the 
charge carriers in graphene nanoribbons (GNRs)^ and 
different gate voltages applied to the two layers of 
bilayer graphene are promising methods of engineering 
band gaps in graphene. 6 Recently, a band gap has been 
discovered in epitaxial graphene derived from graphitiza- 
tion of silicon carbide^ The effect of a substrate-induced 
band gap has also been investigated experimentally 
in epitaxial graphene on nickel^ and theoretically in 
graphene on boron nitride^ 



One of the routes toward tailoring the electronic 
properties of graphene is through the adsorption of 
metals i 10 ^ 1 Alkali metals, in particular, are donors of 
electrons and can be used for the purpose of doping 
graphene to change its carrier concentration. The re- 
lated problem of adsorption of alkali metals on graphite, 
which consists of weakly linked layers of graphene, has 
been studied extensively in the past^ 2 - In a quite recent 
study, a band gap was experimentally observed, and the- 
oretically confirmed, in the surface electronic structure 
of graphite, induced by the adsorption of sodium in a 
(5 x 5) models The similar effect has also been reported 
for the (4x4) models Due to the heavy doping, however, 
the Fermi level is displaced into the graphene conduction 
7T* band so that the band gap lies in the occupied levels. 
The gap- inducing mechanism in these cases involves dop- 
ing as well as interaction between the graphene layers in 
the graphite surface. 

In this paper, we study the gap-inducing effect of 
alkali-metal adsorption on single layer graphene, which 
represents a different problem than graphite surface since 
the gap cannot be due to interlayer interactions. For in- 
stance, first-principles calculations of metal adsorption 
on graphene in the (4 x 4) model 15 do not show a gap 
at the Dirac point in the density of states (DOS). We 
present the case of lithium on graphene compound C@Li 
in detail for several reasons. This is a prototypical sys- 
tem in which the gap opening effect occurs, it is a simple 
system to study theoretically, and it may be accessible 
experimentally. The crucial property here is the struc- 
ture of the system, so that other alkali and alkaline-earth 
metal compounds GqM can exhibit the same effect. For 
instance, a similar band gap has been calculated for cal- 
cium graphite intercalated compound, 16 as well as for 
calcium on graphene, both as CgCaJi Due to doping, 
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lithium or another alkali metal on graphene results in a 
metal and is not directly useful for semiconductor ap- 
plications. It is, however, important to understand the 
mechanism responsible for the gap opening effect. 

An interesting question is which coverages of an alkali 
metal on graphene can induce a band gap in its spectrum. 
Our first-principles density functional theory (DFT) cal- 
culations indicate the absence of a gap in (lxl), (2x2), 
and (4 x 4) models, corresponding to C2Li, CsLi, and 
C32Li, respectively, but its presence in the (\/3x V3)i?30° 
structure, which includes C^Li and CeLi, as well as in 
(3 x 3) coverage CisLi. Explaining the underlying reason 
is the subject of the following sections. 

Our paper is organized as follows. Section HT1 presents 
the results of DFT calculations. Section HTT1 develops the 
tight-binding (TB) model. It includes a description of the 
lattice geometry, from which we clarify the gap mech- 
anism within the TB model. We conclude here that, 
when lithium atoms are arranged above the hollow sites 
of graphene honeycomb lattice in a periodic commensu- 
rate structure to form CeLi, the nearest-neighbor hop- 
ping amplitudes acquire alternating values constructed 
in a Kekule-type structure. This section also includes a 
discussion of the TB Hamiltonian and numerical results 
for the Kekule textured model with two hopping ampli- 
tudes. Section IIVI contains a summary and conclusion. 
The Appendix develops the solution of the TB model. 



II. DFT CALCULATIONS 

Quantitative data for the electronic band struc- 
ture of CeLi were obtained from DFT calculations, 
using a plane wave basis set and pseudopotentials. 
We used the Quantum-ESPRESSO code and norm- 
conserving pseudopotentials,— and performed the cal- 
culations with local-density approximation (LDA) and 
Perdew-Zunger parametrization 19 of the Ceperley- Alder 
correlation functional. 20 We used a lattice parameter 
of ao = 2.46 A for graphene, corresponding to a C-C 
bond length of 1.42 A. Other theoretical calculations^ 
have shown that alkali metals generally prefer to ad- 
sorb on the hollow sites of graphene rather than the 
bridge or top sites. We positioned the lithium atoms 
over the hollow sites of graphene in a (V3 x %/3)i?30° 
commensurate superstructure. We remark, however, 
that lithium has not been shown to coat graphene or 
graphite in this manner J£ As a first approximation and 
in view of the strength of the graphene bonds, initially 
we used the DFT code to only relax the distance of 
lithium atoms from the graphene plane, keeping the co- 
ordinates of carbon atoms fixed. This yielded a distance 
of 1.79 A between the lithium layer and graphene, which 
is a somewhat smaller value, as expected for the LDA, 
than that obtained in other theoretical calculations based 
on generalized-gradient approximation (GGA)i^i. In our 
tight-binding model presented in the Sec. IIIII we use the 
data obtained for this structure. However, since a change 




FIG. 1: (a) Band structure of C6*, i.e., graphene described by 
the unit cell of (y/3 x \^3)R30° structure, (b) Band structure 
of CeLi. A band gap has opened at the neutrality point, 
the degenerate bands along the T-K direction are split, and 
a small band gap has opened at K near the bottom of the 
figure. Horizontal dashed lines indicate the Fermi levels of 
the respective systems. 



in the bond lengths has an important implication for 
opening a band gap, we then used the code to relax the 
positions of the carbon atoms as well as the position of 
lithium atoms. This indeed resulted in a tiny Kekule dis- 
tortion of C-C bond lengths, so that the near bonds to 
Li were contracted by 2.8 x 10~ 4 A, and the far bonds 
to Li were stretched by 5.6 x 1CP 4 A. We then made 
band structure and DOS calculations for CeLi. For re- 
peated images of the systems, we used supercells, with 
c = 12 A for graphene, and c = 15 A for CeLi. A kinetic- 
energy cutoff of 70 Ry was needed for total-energy conver- 
gence. Brillouin-zone (BZ) integrations were made with 
a 6 x 6 x 2 Monkhorst-Pack sampling of /c-space^i with 
Gaussian smearing of 0.05 Ry. For the calculation of the 
DOS we used the tetrahedron method^ 2 - with a finer k- 
point mesh of 36 x 36 x 1 grid. 

The band structures of CeLi, with unrelaxed carbon 
atom positions, and Ce*, i.e., graphene described with 
the unit cell of CeLi, are compared in Fig. [TJ We note 
the fourfold degeneracy of the bands at T, at the Fermi 
level (Dirac or neutrality point), in Fig. [lja), which is 
due to the folding of K and K' points of the graphene 
BZ onto r in the BZ for Ce* structure. We also note the 
linear dispersion of the energy bands in the neighbor- 
hood of the Dirac point. In Fig. [IJb), we see a number 
of changes that appear in the band structure of CeLi. Of 
interest for our study is the gap opening of E g = 0.39 eV 
at the neutrality point. For the CeLi with fully relaxed 
positions, including those of carbon atoms, the gap open- 
ing is Eg = 0.41 eV, and otherwise the band structure 
is quite similar to the unrelaxed case. This means that 
the change in C-C bond lengths is responsible for only 
about 5% of the band gap. We also note the lifting 
of other degeneracies, at K and along the T-K direc- 
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tion. Other changes that can be readily noted, as a re- 
sult of lithium adsorption, include the appearance of the 
parabolic-shaped band of lithium, with its minimum lo- 
cated at F, at 0.92 eV above the neutrality point. The 
lithium band weakly hybridizes with the graphcne con- 
duction bands, resulting in some shifts and kinks in the 
bands of graphene. There is a substantial charge trans- 
fer of > 0.2e per Li atom from lithium layer to graphene, 
with every carbon atom receiving the same charge of one- 
sixth of this value, i.e., > 0.03e. As a result, the Fermi 
level is raised by 1.5 eV relative to the neutrality point of 
graphene. A rigid band model is clearly not a complete 
description, but it can be used to describe the charge 
transfer, which is another important aspect of this prob- 
lem. 



III. TIGHT-BINDING MODEL 

A. Lattice geometry 

The band gap of Cghi is a consequence of the special re- 
lationship between its structure and that of the graphene 
substrate. The Bravais lattice of CeLi is hexagonal, as 
it is for graphene, but with a times larger lattice pa- 
rameter and rotated by 30°, as shown in Fig. [2] It is 
convenient to describe the base vectors of the hexago- 
nal lattice in a standard way, with the lattice parameter 
used as unit of length, so for CgLi we use the coordinate 
system with its cc-axis taken along ai in Fig. [21 




FIG. 2: (Color online) Lattice structure of lithium on 
graphene CeLi. Carbon atoms on sublattices A and B are 
depicted as filled and open black circles. Lithium atoms, in 
magenta, occupy positions above the hollow sites. The dashed 
hexagon is the Wigner-Seitz primitive unit cell of CeLi, con- 
taining six numbered C atoms and one Li atom. Two types 
of bonds, two-thirds red and one-third blue, are distinguished 
due to the presence of Li atoms. On the lower left, the rhom- 
bus defined by ci and C2 and enclosing two carbon atoms is 
the primitive unit cell of graphene. The 30° rotated rhombus, 
defined by ai and a.2 and having y/3 times larger side, is the 
unit cell of CeLi. In the upper left part are shown the vectors 
pointing from a carbon atom on sublattice A to its three near- 
est neighbors on sublattice B. The inset shows the reciprocal 
space. The vectors bi,2 and di,a are sets of reciprocal lattice 
vectors of CeLi and graphene, respectively. The hexagons in 
the inset are the Brillouin zones. The region defined by YKM 
is the irreducible wedge of graphene Brillouin zone. 



ai = (l,0)a, 



a 2 = - 




(1) 



where a = aoV3. [Thus the same base vectors [Eq. (J])] 
describe graphene, i.e., Ci and c 2 in Fig. [21 with a = ao 
and with the x axis taken along ci.] The corresponding 
reciprocal-lattice vectors of the hexagonal lattice are 



1 

V3 



2 \ 2tt 



We denote the reciprocal-lattice vectors of graphene, 
which are also given by Eq. ^ relative to its own sys- 
tem, by di^- The two sets of vectors hi 2 and di. 2 are 
related by a 30° rotation and a v3 change in scale, i.e., 
|b| = |d|A/3. 

As shown in Fig. [21 the K points of graphene are lo- 
cated at 



K: idx + ids^bx, 



K': -id 1 + ^d 2 ^b 2 , (3) 



i.e., the K points of graphene coincide with reciprocal- 
lattice points of C6*- It also follows from Eq. (J3J) that 
the K points of graphene belong to the reciprocal lattice 
of the (3 x 3) structure. Evidently, (n\/3 x nV3)R30° 
and (3n x 3n) constructions, with n a positive integer, 
will also share this property. The existence of this kind 



of relationship is the basic reason why a gap opens in 
C X M if a; is a multiple of 3, since in these structures the 
K points of the underlying graphene become coupled and 
the perturbation caused by the foreign atoms removes the 
degeneracy. The wave functions with wave vectors corre- 
sponding to K and K' of graphene mix to form different 
standing waves. One standing wave piles up electronic 
charge in hexagonal cells occupied by Li, and depletes 
it from the hexagonal cells devoid of Li, and the other 
standing wave does the opposite. (The commensurate 
structure of C@Li breaks a Z% symmetry.) These stand- 
ing wave states then experience different potentials which 
qualitatively explains the origin of the gap. 

Each carbon atom of sublattice A is connected to three 
carbon atoms of sublattice B by the r vectors, as shown 
in Fig. [21 These vectors in the coordinate system of CeLi, 
needed in our tight-binding calculations, are given by 



Tl,3 



1 1 



T-2 = 



(4) 



We now examine the system from the viewpoint of the 
tight-binding model. The presence of the lithium atoms, 
in CeLi, above the hollow hexagonal sites can be imag- 
ined to have two effects, within a simple tight-binding 
model: (a) a contribution (Hartree potential) to the on- 
site energies of the C atoms, and (b) a possible change in 
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the hopping amplitudes depicted as C-C bonds in Fig. [2] 
From an inspection of the Wigner-Seitz cell in Fig. [21 the 
C atoms are seen to occupy equivalent positions with re- 
spect to the metal atoms since each is at a vertex shared 
by one metal-filled and two empty hexagons of the hon- 
eycomb lattice. Therefore all C atoms receive the same 
charge transfer from Li atoms, as we have also seen in our 
DFT calculations, and the symmetry breaking cannot be 
attributed to on-site energies. This reduces the number 
of parameters and makes a TB model description much 
simpler for CgLi than for more dispersed compounds. 
The bonds, however, occupy two kinds of positions in 
a Kekule construction, two-thirds of them (colored red) 
between a filled and an empty hexagon, and one-third 
(colored blue) between two empty hexagons. We con- 
clude that there are two different hopping amplitudes 
corresponding to the two kinds of bond positions. Our 
DFT calculations showed that changes in bond-length 
account for ~ 5% of the energy gap. The modulation in 
hopping amplitudes, which are matrix elements of the 
Hamiltonian between carbon 2p z orbitals, is therefore 
caused mainly by different potential energies and elec- 
tron concentrations in the regions of red and blue bonds, 
with a small contribution due to bond-length distortion. 
Of course, both variations in potential energy and the 
distortion of bond lengths are caused by the presence of 
lithium ions, and their effects on the band gap add with 
each other. 



B. TB Hamiltonian 

The 7r bands of graphene are described well by a Hamil- 
tonian of spinless electrons hopping on the honeycomb 
lattice of graphene, 



H 



E eiC * c * 



c]ci. 



(5) 



Here a (cj) annihilates (creates) an electron at site i, €i 
are on-site energies, which are all equal and therefore can 
be set to zero, and Uj are the hopping amplitudes. 

The nearest-neighbor TB model of graphene has a gap- 
less band structure, but a Kekule modulation of the hop- 
ping amplitudes couples the K and K' points, leading 
to a mixing of degenerate states and opening of an en- 
ergy band gap. 23 Following our discussion above, we al- 
low the hopping amplitudes tij to have two values, t\ 
for two-thirds of the bonds on the hexagonal cells be- 
neath a metal atom, and ti for the remaining one-third 
of the bonds. Here we are interested in describing the 
gap, and therefore we exclude the metal atoms from our 
TB model, except for their effect on modulating the hop- 
ping amplitudes. They can be included in a more general 
treatment j2a but this is not needed for describing the gap 
of C 6 Li. 

For clean graphene t\ = t% = t. There are two 7r 
bands corresponding to the two-atom basis, and the band 
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FIG. 3: Band structures of graphene in the Brillouin zones of 
(a) the (1 x 1) structure, and (b) the (%/3x V3)R30° structure. 



energy can be written as 



E±(k) = ±t 



i=l 



(0) 



where the vectors Xj arc defined in Fig. [2] When we 
use the unit cell of the (s/3 x \/3)R30° superlattice to 
describe graphene, we will have six w bands since there 
are six carbon atoms within this unit cell. The addi- 
tional bands are derived from Eq. ([6]) as E±(k — bi) and 
E±(k — \)2), with the vectors b 1]2 given by Eq. ©, which 
result in the folding of K and K' points onto the T point 
of BZ, because of relation ([3]). The energy bands from 
these analytic formulas are plotted in Fig. [3j ft is seen 
that some degeneracies are present at the new K points 
[Fig. [31(b)] as well as along TK and MK lines. The linear 
dispersion near the Dirac point is shown in the Appendix 
to be E±(k) = ±Hvf\^\, where, vp = \^3aQt/2h is the 
Fermi velocity, t is the nearest-neighbor TB hopping pa- 
rameter, and ao is the graphene lattice parameter. 



C. Numerical results 

We must verify that the band gap of CeLi can be de- 
scribed by a textured TB model with two different hop- 
ping amplitudes, t\ ^ t 2 . In this case, we obtain the band 
structure by numerical diagonalization of the Hamilto- 
nian [Eq. ([5])] in momentum space. The method is de- 
scribed in the Appendix, where we also derive an analytic 
expression for the energy dispersion in the neighborhood 
of the gap as 



E± = ±l h'v F 2 k 2 



where 



E„ = 2\h 




(7) 



(8) 
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for the textured TB model. The Fermi velocity is given 
in terms of the TB parameter as before, except that now 
an average value, 



must be used. 

From the slope of the linear dispersion near the neu- 
trality point of Fig. [Ija) , we can find the Fermi velocity 
and then the TB parameter t = 2.72 eV, in good agree- 
ment with the commonly used range for this parameter^ 
From Fig. [2](b) , we can find the renormalizcd Fermi ve- 
locity vf and thus the hopping amplitude t by fitting 
Eq. ([7]) to the band-structure data in the vicinity of the 
gap. We found that the value of t is slightly larger for the 
upper branch than that for the lower branch, by ~ 0.1 eV, 
but that the change in the average value is much smaller 
than this amount. In the following, we neglect the renor- 
malization of the Fermi velocity and the average hopping 
amplitude. 

Given the values of t = 2.72 eV and E g = 0.39 eV, from 
first-principles calculations, we can determine t\ and 
from Eqs. © and ©. Thus we find t\ = 2.79 eV and t 2 = 
2.59 eV. Interchanging t\ and t% is also a solution since 
Eg depends on the absolute value of their difference. To 
distinguish between the two, we plot the band structures 
corresponding to both results in Fig. 0J The band gaps 
are identical in Figs. BJa) and|3Jb), but a characteristic 
difference can be seen in the way the threefold degeneracy 
is lifted at K. Comparing with Fig. [TJb), it is seen that 
Fig. Ufa) has the same qualitative feature, i.e., the same 
shape of band gap, and the lifting of degeneracies along 
YK and at the K point. This shows the validity of the 
Kckulc textured model for describing lithium adsorbed 
graphene, as well as showing that for CgLi, t\ > t%. In 
fact, Fig. Hlb) corresponds to our DFT calculations of 
the band structure of CaLi (not shown), which has the 
same lattice structure as CeLi, and is its complement. 

A useful quantity that can be calculated from energy 
bands obtained is the density of states per unit cell given 
by 

9 (£) = 2aW 7^8{E-E n (k)}, (10) 

n=l"' BZ * 77 ' 

where A is the area of the Cq* unit cell. For the TB DOS 
we used a 60 x 60 /c-point mesh and the triangle method, 26 
which is the two-dimensional version of the tetrahedron 
methodS 2 - that we also used for the ab initio calculations. 
The density of states calculated by first-principles DFT 
and our TB model compare favorably in the vicinity of 
the gap, as shown in Fig. The extra densities of states 
in the energy range above the gap in Fig. [5ja) are due 
to the band of the lithium atoms, which are not included 
in our tight-binding model. The other main difference 
between the two results is the positions of the Van Hove 
singularities which are closer to the gap in first principles 
calculations. 



(a) (b) 
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FIG. 4: Band structures of graphene from the TB model with 
two hopping amplitudes: (a) t\ = 2.79 eV, ti = 2.59 eV, 
and (b) ti = 2.59 eV, ta = 2.79 eV. Note that the three-fold 
degeneracies at K, e.g., the lower left one, are lifted so that in 
(a) we have first a non-degenerate level and above it a doubly 
degenerate level, and in (b) we have just the opposite order 
of the degenerate and nondegenerate levels. 



IV. SUMMARY 



We have studied the opening of a band gap induced 
in single layer graphene by the adsorption of a sub- 
monolayer of an alkali metal. A band gap induced by 
the adsorption of alkali metals was previously observed 
in graphite, and its origin was explained to be due to 
charge transfer to substrate and interlayer interaction of 
graphene layers. In single layer graphene, where inter- 
layer interaction does not exist, the origin of an alkali- 
metal induced band gap is in the coupling of K and K' 
points which is only possible for appropriate superstruc- 
tures. The required property exists in the superstruc- 
tures (3nx3n) and (nv / 3xnv / 3)i?30 o , with n as apositive 
integer. For other coverages on graphene we expect the 
bands to remain degenerate at the K points of graphene, 
resulting in a gapless spectrum. A tight-binding model 
with alternating hopping amplitudes forming the Kckule 
construction was shown to describe the characteristic fea- 
tures of the band structure of CeLi. From our DFT calcu- 
lation, we concluded that the modulation of the hopping 
amplitudes is mainly due to variations in charge density 
in the graphene plane due to the attraction caused by 
the Li ions. A tiny distortion of bond lengths additively 
contributes ~ 5% to the band gap and, thus, enhances 
the modulation in hopping amplitudes. In conclusion, 
we showed that the band gap opening in the lithium on 
graphene compound CgLi is consistent with the nearest- 
neighbor tight-binding model of graphene with Kekule 
modulated hopping amplitudes. 
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E(eV) 

FIG. 5: (Color online) (a) DOS of C 6 Li calculated by DFT 
code. The vertical red line indicates the position of the Fermi 
level, (b) The corresponding DOS for graphene calculated 
within the tight-binding model. The position of Van Hove 
singularities are displaced relative to the DOS of (a) and the 
extra densities corresponding to Li atoms are not included in 
(b). 
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APPENDIX: SOLUTION OF TB MODEL 

In this appendix, we summarize our method of solv- 
ing the Kekule textured nearest-neighbor tight-binding 
model [Eq.[5]. First, we expand the eigenfunctions of our 
general periodic system, characterized by crystal momen- 
tum k, in terms of the 2p z orbitals, \nj), localized on the 
jth atom in the nth unit cell, as 

N M 
ViV n=lj=l 

Here N is the number of unit cells, and M is the number 
of carbon atoms in the basis, r n j are position vectors of 
the atoms, and the amplitudes Aj form the components 
of the eigenvectors to be determined by diagonalizing the 
Hamiltonian. For graphene M — 2, and the Hamiltonian 
matrix in momentum space is given by 

flW =-'k«-" • El fl' (A - 2) 

which is readily diagonalized to give Eq. ©. 

For the textured model, M = 6. Referring to Fig. [2] 
for numbering of the basis atoms, and the t, vectors, 
we can write the 6x6 Hamiltonian matrix. For exam- 
ple, to obtain the first row we observe that atom 1 has 
nearest neighbors 2, 4, and 6 (due to periodicity), and is 
connected to them by T3, t 2 , and Ti, respectively, with 
hopping amplitudes t\, t\, and t 2 . Thus we find 



ff(k) = - 










tie 1 ^ 







tie lW ' T2 







i 2 e ik - Tl \ 


tie 


-ik 


T3 





tie 


— ik*Ti 





t 2 e 


-ik-T 2 












tie lk - Tl 







t 2 e lk - T3 







tie* -7 " 2 


tie 


-ik 


V2 





t 2 e 


-ik--r 3 





tie 














t 2 e lk - T2 







tie lk - Tl 







tie lk - T3 


\t 2 e 


-ik 


TI 





tie 


-ik'-r 2 





tie 


-ik-T 3 


J 



(A.3) 



For general k € BZ, we diagonalize Eq. (|A.3[) numerically, but we obtain an analytic expression for the neighborhood 
of the gap using perturbation theory. 

At k = 0, the unperturbed Hamiltonian, case of t\=t<i = t, has four of its eigenvalues equal to 0, with the other 
two being ±3t. The corresponding eigenvectors are easily written down by noting that they are related to the K, K' , 
and r points of graphene described by its primitive unit cell, 



1 








1 

7! 



id 

UJ* 





1 

7! 





UJ 





1 

7! 





u* 


LJ 



1 

7! 



±1 
1 

±1 

V±i/ 



(A.4) 
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where u = e l27r / 3 . 

Next we calculate the Hamiltonian matrix in the fourfold degenerate subspace of the zero eigenvalue, spanned by 
the first four vectors of Eq. (|A.4p , and expand the result near k = to obtain (cf . Ref . 23) 



W(k) = 



/ hvF(k x — iky) t"i—t\ 
hvp{k x + iky) t2~ti 

t<x — 1\ —hvF(k x — iky) I ' 

\ ti — 1\ —hvp(k x + ik y ) 



(A.5) 



where vp — y3aot/2H as for unperturbed graphene, but with t — (2ti + t<z)/Z. Diagonalizing Eq. (|A.5|) . we find 

Eq. Q. 
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